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Abstract 

Stability of homogeneous precession of spin due to parametric ex- 
citation of spin waves is considered as the explanation of the "catas- 
trophic relaxation", that is observed in the superfluid "^He — B . It is 
shown, that at sufficiently low temperatures homogeneous precession 
of spin becomes unstable (Suhl instability). At zero temperature in- 
crements of growth for all spin wave modes are found. Estimation of 
the temperature of transition to the unstable state is made. 

1. The use of pulsed NMR is based on the investigation of homogeneous 
precession of spin in a constant magnetic field. Spin precession induces the 
free induction signal (FIS), which is registered in the induction coil. Preces- 
sion of spin in the superfluid B-phase of ^ He has its specifics. At temperatures 
T > 0,4Tc, where Tc is the temperature of superfluid transition, FIS exists 
anomalously long - many times longer than the time of dephasing of spin due 
to a residual inhomogeneity of d.c. field. At T ~ 0,4Tc occurs transition 
into the other regime, when on the contrary FIS disappears quickly. This 
fast decay of precession was first observed in ref. P and since that is re- 
ferred as catastrophic relaxation. While anomalously long FIS was explained 
a long time ago, even a qualitative explanation for the catastrophic relax- 
ation is missing. The decay of homogeneous precession was demonstrated by 
numeric simulation of equations of spin dynamics p!. The simulation was 
made in the restricted geometry and the authors of simulation attribute the 
main role in the destruction of precession to the walls, i.e. they consider the 
mechanism of destruction to be surface. 
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In the present paper the explanation of catastrophic relaxation is sug- 
gested, which is based on a bulk effect. It is the instability of homogeneous 
precession with respect to decay into parametrically excited spin waves with 
opposite wave vectors (Suhl instabihty [Sj). 

Fast decay of precession was also observed in the uudd -phase of solid 
^ He and it was explained by the onset of Suhl instability pi . However, the 
quantitative interpretation of the results, concerning instability of precession 
in the cited work is based on the modification of theory Ref. 01, developed for 
the continuous NMR and which is therefore applicable only for small tipping 
angles. In our analysis this restriction is not used and it can be applied for 
the arbitrary angles between spin and magnetic field. 

2. In order to describe motion of spin we use expression of Hamiltonian of 
Leggett, which is written in terms of Euler angles a , jS , $ = + 7 (z-axis is 
oriented opposite to the direction of d.c. magnetic field Ho ) as coordinates 
and canonically conjugated momenta P = Sz — S(; , , 5^ , where 5*2 — 
is projection of spin onto z-axis, Sq — its projection onto axis ( = Rz and 
5*^ — is projection on the line of nodes (see for details [6 ). We choose units 
of measurement so that x = (7 = 1 , where x — is magnetic susceptibility 
per unit volume of ^He — B , and g — is the gyromagnetic ratio for nuclei 
of ^He ; after that spin has dimensionality of frequency and energy of the 
frequency squared. Using variable and units mentioned above one can write 
the Hamiltonian in the form: 

(1) 

where ul - is Larmor frequency corresponding to the d.c. magnetic field, 
Fy - gradient energy, $) - dipole energy. Ud for ^He — B depends 

only on two variables (3 and $ , that justifies the choice of $ as a variable, 
when the dipole energy is essential. Gradient energy for ^He — B can be 
written as: 

-^V = "^[cl^ik^ir, - (Cj| - c\){5i^5krj + 5ir^5k^)]uJi^UJkn, (2) 

where 

oji^ = — sin /? cos 7 + sin 7, 

u;2^ = sin/5sin7 -|- C0S7, (3) 
= a,^cos/3 + 7,^, 

= ^ etc.,c| c\ — are squared velocities of two types of spin waves 
("longitudinal" and "transverse"). In what follows we choose units in such 
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a way as to Cy = 1 , then wave vectors entering the equation will also have 
dimensionality of frequency. Equations of motion, that are generated by 
Hamiltonian have the form: 

da _dH dP _ dH d dH 
dt dP ' dt da dx^ da^^ ' 

etc. 

System of equation (j3)) has spatially uniform stationary solution describ- 
ing precession of spin in the stationary magnetic field at Q < (3 < 9q = 
arccos(— |) : 

a = uJit + ao, 7 = -uJit + $° - «o, 
P(°) = cjl(cos/? - 1), Sf = 0, = ulcosP, (5) 

cos $(o) = (l-cos/3W)/(l + cos/3^°^). 

It is convenient to introduce ■?/' = « + uj^t instead of a and at the same time 
to transform the Hamiltonian H = H + ojlP , so that ^ = . 

Let us now obtain equations for the small deviations from the stationary 
solution: 

5V'(r, t) = V - V'^"^ (6) 

etc. 

In zeroth approximation on the small deviations the gradient energy has 
three groups of terms: "stationary" - with time- independent coefficients, and 
two "oscillating", corresponding to Larmor and doubled Larmor frequencies. 
Without the loss of generality we consider perturbations propagating in yz 
plane: 



Fsjst = ^[^i^lil - /isin^/jW] + (1 _ ^)5/52 ^ ^^2^ ^ 25^,,57,,cos/3(°) + 

- 2/icos2/5(°)] + 5pl + 57^(1 - 2/i) + 25^,,57,,cos/5(°)(l - 2/i)], (7) 

Fvosci = -l^[S'ip,yS'ip,z sin 2/5^°^ sincuit + 
+ {Silj^yS'-f^z + Sip^zSj^y) sin/5*-°'' sint^Lt + [Stp^yS/S^z + ^V'.^'^As/) cos/5''°^ costUit + 

+ {6 P 6-f ^z + SP,zS7,y) cos ULt], (8) 

F^osc2 = -| [sin' P^^'^S^/jI cos 2iULt - 6/31 cos 2c^it + 

-|- sin 
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Oscillating terms are proportional to /i = 1 — c^/c|. We assume that fi 
is small, and consider oscillating terms in the equations of motion as small 
perturbations. Actually fi is not very small ~ 1/4 in a vicinity of 
Tc), however, this approximation gives satisfactory results. More precise 
criteria of application of such approximation will be formulated in a process 
of solution. 

In order to use the theory of perturbation we write linearized system of 
equations of motion in the form: 



Hi 



Mo + Vm X, (10) 



where 



X 



( 5i) \ 
6(3 
5$ 
6S, 
SSp 

V J 

Matrix operator Mq includes all time-independent terms, and oscillating 
terms are collected in matrix operator V{t) , which is proportional to /i and 
therefore is considered as a perturbation. Equation of zero order approxima- 
tion on perturbation: 

f^A^oX, (11, 

gives us dispersion laws for the three branches of spin waves: 



(12) 



lu;l + e+'^iJL\/u;l + 4k^ 



and eigenvectors, corresponding to each oscillation branch Xj(/c), i = 1, 2, 3 . 
Recall that we are speaking about the spin waves propagating against the 
background of homogeneous precession, so the mentioned oscillations of spin 
are different from the usual spin waves, which are created by small deviation 
from equilibrium orientation. It is assumed here that cu^ 3> where Q is 
the frequency of longitudinal oscillations. 
Solution of the system (fTT|) is given by: 



X(t, k, y, z) = ^{ai_Xj(fc) exp{—iLJi{k)t + ikr) 



2=1 



aj+X*(A;) exp{iuji{k)t + zkr) 
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a*_X*(fc) exp{iLJi{k)t — ikr) + 
a*^^i{k) exp{—iuji{k)t — ikr)}, (13) 

where Oj - are constant coefficients, k = {ky, kz) and r = (y, z) are 2D 
vectors. When V{t) is taken into account ^{t, k,y, z) , given by equation 
(fT^. does not satisfy the equation of motion (fTUj) . The ffist order approx- 
imation on V can be obtained by the method of averaging of the classical 
mechanics. Substituting zero order approximation (fT!?jl into formulas (jHI),© 
we can see that terms of the first order on fi will not vanish after time aver- 
aging only if there are resonance relations between ul and eigenfrequencies 
uji{k) : Wi{k) = ojl and uJi{k) = ujl/2. As it is seen from equations (12) 
for all branches of spin waves there exist k which are satisfies such reso- 
nance conditions. Vicinities of these wave vectors are "dangerous" for the 
appearance of instability. From the same formulas it follows that for the 
different branches resonance conditions are satisfied with different values of 
wave vectors, and therefore each branch can be considered independently. 
To find solution nearby the resonance frequencies we will use the standard 
procedure, when solution in the main approximation is sought in a form: 

X.^^\t,k' ,y,z) = a,„(t)Xi(A;')exp(-iu;2^t + «k'r) + 
ai+{t)X*{k') exp(iu;2^t + ik'r) + 
a*_{t)X.*{k') expiiJ^^h - zk'r) + 
a*^{t)Xi{k') exp{-iJJ^h - ik'r), (14) 

/ = 1,2, where w]^ , - is one of the resonance frequencies (cujj^ = 

(2) ' 

c<j)j. = cul , we will suppress index / in the nearest formulas for brevity), k - is 
wave vector nearby resonance frequency for the i-mode, aj(t) - are "slowly" 
varying functions of time, i.e. d± ~ fia± . Terms, that have frequencies 
differing from on integer multiple of 2uJi — 3ljji,5ujji,7u!r appear in 
the next orders on /x . 

Substitution of solution (fT^ into system (fTUjl gives us: 

X„ + X^^ = M(A;')X(t, k\y, z) + V(k', t)X(t, A;', z), (15) 

where 

Xq = dj_(t)Xj(A; ) exp(— icj/jt + ik r) + 
aj4.(t)X*(/c ) exp(zco'/jt + ik r) + 
d*_(t)X*(fc') expiiujRt - ik'r) + 
d*^(t)Xi(A; ) exp(— icij/jt — zk r), (16) 
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X,^^ = -iujjiai^{t)Xi{k') exp{-iujRt + ik'r) + 
iujRai+{t)X*{k ) exp{iujjit + ikr) + 
iuo iia*_{t)'X*{k ) exp{iu!jit — ikr) + 
—iu!Ra*_^_{t)Xi{k ) exp{—iujRt — ik r). (17) 

Taking into account that Xj(fc') and X*(A;') are eigenvectors of M , cor- 
responding to the frequencies uji{k ) and —uJi{k ) one can rewrite equation 
(US} as: 

X, + X^^ - X^, = V^(k', t)X(t, fc', I/, z), (18) 

or 

X, - ^X^^ = f (k',t)X(t, A;',!/,^), (19) 

where e{k') = uj{k!) — lur . 

Let us multiple the last equation by exp(— zk r) and take integral over 
volume. As the result a*_ and a*^ vanish. Expressing cosines and sines in 
terms of exponents: 

X« - ^Xf), = ivl'\^) exp(2.c.«t) + V^'\l^) exp(-2.u;«t) + 

+Vi')(k') exp(2^u;g)t) + vi')(k') exp{-2tJ^h))X^'\t), (20) 

one obtains the sum of terms with different powers of exponents. Since 
V{t) contains cosines and sines of 2u!^J^t , coefficients ai-{t) and cij+(t) are 
related by exponents with the same powers ±iuj^^t . The resulting equations 
are multiplied by exp{±iujfj^t) and averaged over rapid oscillations. Finally, 
after making projection of equations on eigenvector of i-mode one obtains 
system of two differential equations of first order which relates aj_(t) and 
ai+{t) : 

d,+ (t) - z£«(A:')a,+ (t) = ^ Vi2 ' ^ --(^)^ (21) 
dUt)+ie^'\k')a,4t) = <X.^£(k)X.> ^^^^^^_ ^22) 
System ()21|) has solution proportional to exp (A^'^t) , where A^'^ is defined 

by: 

,(0 , ±i (< x.^>;"(k-)x. >< x;v:!"(k-)x- > _ ^ ^^^^ 

2 \ |Xj I / 
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Resonance corresponds to the value of /c' , when e^^\k') = . In a region 
of k close to resonance expression in the brackets is positive. Then one of 
the values of A^'-* corresponds to the growth of amplitude of oscillations, i.e. 
development of instability begins. 

3. Let us consider all possible cases of resonances. For each mode we will 
write: law of dispersion, eigenvector of this oscillation and increment, which 
is obtained on the condition of resonance. 

First mode. Law of dispersion: 



(24) 



Eigenvector: 



Xi.{k) 



Resonance at the frequency 





1 

—iLJi{k) cos j3 


V iuji{k){l - COS/9) J 



2 



± 



Increment of growth: 



(1) 



/i- 



ul sin — 2 cos /3| 
T ' 2 cos2 /5 - 2 cos /5 + 5 



■ sin 25, 



(25) 



(26) 



where 5 - is the angle between direction of the wave vector and z-axis. 
Maximum increment corresponds to the direction: 



= ^. (27) 
increment vanishes in the case of wave vector directed 



As it is seen from 
along y-axis. 

Resonance at the frequency ■ 

k = ±ujl- 

In zeroth order approximation on dipole frequency we have: 

Af ) = 0. 



:28) 



(29) 
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Finite increment appears when the dipole terms are taken into account in 
the equations of motion. In this case resonance condition is satisfied by the 
wave vector: ^ 

k' ^UL-^il + ^cosP)—, (30) 

and eigenvector has corrections of the order of Vt^ /uj\ . With these corrections 
the increment will be equal to: 



'(2) _ /i . 2 ^(1 +4COS/3) 



1/2 



oujl (1 + cosp)-^/^ 



|l-2cos/3|; 



sin^ 5 



2 + 2 cos2 13 -2 cos 13 
Maximum increment corresponds to the direction: 



'(0) _ TT 



(31) 



(32) 



Second mode. Law of dispersion: 

o 1 



Eigenvector: 



i sin P 



1 



V 



1 — cos P 

U2[k) 

zsm^ • ( ^^) 

(sin^/?) 



<^2(A;) 



Resonance at the frequency ^ : 



A; 



± 



"\/3a;L 



(33) 



(34) 



Increment of growth (to zero order approximation on dipole energy) 

A« = 0. 

Resonance at the frequency ujl (taking into account dipole energy): 

2 Vt^ 
k' ^V2ul- —{l-cosP) — . 

15 ul 



(35) 



(36) 
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Increment: 



^[i) ^ _p — sin^/3(i-cos/?)|l-4sin2/3|- 



5 lul '■^'12-2cos/3-17cos2/3 + 8cos4/3' 

(37) 

has its maximum for ^2 = vr/2. 
Third mode. Law of dispersion: 



^3 = l^l + k' + \ujL^ul + Ak\ (38) 

Resonance at the frequency ^ is not possible because the frequency of 
this mode is larger than ujl for all k . Resonance at the frequency uo^ is 
also not possible because near /c = A is imaginary. 

As it is seen from the Fig. there exists positive increment at least for 
one of the modes of oscillations for all tipping angles. The main role plays 
maximum increment that is found for the first mode in the case of resonance 
frequency ujl/2 . 

4. The obtained results are correct at T = . At a finite temperature spin 
waves damp. This leads to appearance of temperature threshold of instability. 
We can take into account small damping by substituting new complex law of 
dispersion into formula for increment ()23|1 . Here we should replace e^^\k Y 
by \e^^\k )p. To estimate temperature threshold of instability one should 
use the law of dispersion with the corrections for damping: 

cu;2(A;) = ujl{k) - 2iD{T)ujLk'' + 0(A;"), (39) 

where D{T) - is coefficient of diffusion. We substitute into corrected 
formula (j^Hj) with , which satisfies the resonance condition Re{uj{k)) = 
. After substitution and with account that Re{uj{k)) = uj^^ we arrive at: 

>^'HT) = ^{xL.iT = 0)-D\T)'^^f^ . (40) 

This formula determines temperature below which instability sets on: 

\(i) ,,(0 

DiT) = 7;. (41) 

We estimate coefficient of diffusion at T = 0.4rc by using (jH]) ( Tc - is the 
temperature of superfluid transition). For the first mode of oscillations and 
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for the resonance frequency ujl/2 at tipping angle 90° and for the pressure 
20 bar: 

2A(^) cl 

D{T) = "7 " = 0.027 W/s, (42) 

C|| = 2050^^1 - T/Tc ^ 1600 sm/s, c = 3/4 , ujl = 2.9 ■ 10^ s . This result 
can be compared with the experimental data for the transverse coefficient of 
diffusion in ^He — B [8\. In the cited work the value of transverse coefficient 
of diffusion is approximately equal to 0.03 sm^/ s for the pressure 20 bar 
and for the Larmor frequency 2.9 ■ 10^ s~^ . Thus the estimated "critical" 
coefficient of diffusion is close to the measured one at 0.4Tc . 

5. It follows from the given analysis that at sufficiently low tempera- 
tures homogeneous precession of spin in ^ He-B is unstable because of Suhl 
mechanism. Interaction between precession and spin waves appears mainly 
because of the anisotropy of spin wave velocities. Estimation of decay time of 
precession as inverse maximum increment of growth of spin wave amplitude 
gives result which does not contradict to the measured value at the lowest 
temperatures. Estimation of the temperature threshold of the onset of insta- 
bility using the available data about the value of spin waves damping falls 
into the temperature interval in which transition from stationary precession 
to catastrophic relaxation is observed. This allows to consider Suhl instabil- 
ity as the probable reason of the observed catastrophic relaxation. In order 
to make proposed here explanation of catastrophic relaxation quantitative 
one should describe more precisely spin waves damping taking into account 
direction of propagation. This work is in a progress. 

This work is partly supported by RFBR (grant 04-02-16417) Ministry of 
Science and Education of the Russian Federation and CRDF (grant RUPl- 
2632-MO04). 
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Figure 1: Dependence of increments for two spin waves modes: (1) - 

AS^)(/3W,5r)/^L (I2SD,(2)- A?)(/3W,5r)/^L (EU, (3) - A?)(/3W, ^^"V^l 
(jnZI), on tipping angle /5o in interval < /3 < = arccos(— , VLjuji^ = 1/2 
for (2) and (3). 
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